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Abstract. We adapt the bialgebra and Hopf relations to expose internal 
structure in the ground state of a Hamiltonian with Z2 topological order. 
Its tensor network description allows for exact contraction through simple 
diagrammatic rewrite rules. The contraction property does not depend on specifics 
such as geometry, but rather originates from the non-trivial algebraic properties 
of the constituent tensors. We then generalise the resulting tensor network from 
a spin-i lattice to a class of exactly contractible states on spin-S degrees of 
freedom, yielding the most efficient tensor network description of finite Abclian 
lattice gauge theories. We gain a new perspective on these states as examples 
of two-dimensional quantum states with algebraically contractible tensor network 
representations. The introduction of local perturbations to the network is shown 
to reduce the von Neumann entropy of string-like regions, creating an unentangled 
sub-system within the bulk in a certain limit. We also show how perturbations 
induce finite-range correlations in this system. This class of tensor networks is 
readily translated onto any lattice, and we differentiate between the physical 
consequences of bipartite and non-bipartite lattices on the properties of the 
corresponding quantum states. We explicitly show this on the hexagonal, square, 
kagomc and triangular lattices. 
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1. Introduction 

The inherent difficulty in describing quantum many-body systems is a consequence of 
the "curse of dimensionality" — the cost of describing a state grows exponentially with 
the number of degrees of freedom. Progress has been made on this front, exploiting 
results from quantum information theory which constrain physically relevant states 
to a small corner of Hilbert space. In one dimension the problem has been effectively 
solved for short-range interactions with matrix product states (MPS) [1] and the 
density matrix renormalisation group (DMRG) [2, 3, 4, 5], but in higher dimensions 
DMRG becomes computationally unfeasible. To tackle this, projected entangled pair 
states (PEPS)/tensor network states (TNS) [6, 7, 8, 9] directly generalise MPS, while 
tree tensor networks (TTN) [10] and the multi-scale entanglement renormalisation 
ansatz (MERA) [11, 12, 13, 14] instead make direct use of a hierarchical structured 
tensor network. 

Tensor networks are capable of describing quantum states ranging from those 
with short-range entanglement (such as symmetry-breaking states) up to topological 
phases of matter. However, while the tensor network approach is highly versatile in 
terms of the quantum states it can describe, its utility is threatened by the difficulty 
of contracting the tensor network to obtain reduced density matrices and expectation 
values. Approximate methods exist to perform this contraction. PEPS were first 
introduced along with an algorithm based on interpreting the network contraction in 
terms of MPS and matrix product operators (MPOs) [15]. The tensor renormalisation 
group (TRG) [16, 17] method has also been presented as a means for performing this 
contraction approximately in polynomial time, and has been demonstrated to provide 
an accuracy beginning to approach that of quantum Monte Carlo (QMC) methods 
where applicable [18, 19]. Other methods to efficiently contract tensor networks 
include corner transfer matrix methods [20] and a coarse-graining method explicitly 
including and exploiting a local symmetry for computational gain [21]. 

By contrast, our interest here is in developing a class of quantum states having 
TNS descriptions which are exactly contractible through simple diagrammatic rewrite 
rules [22, 23]. Already it is known that MERA forms an efficiently-contractible class of 
PEPS [24], and that certain states (such as the MERA representation of the toric code) 
arc fixed points under an entanglement renormalisation procedure [25, 26]. However, 
the simplicity of our TNS contraction docs not rely on the underlying network 
topology, on tensors constrained to be unitaries or isometries, or on being a fixed point 
of an abstract renormalisation procedure. Rather, these network rewrite properties are 
graphical rules which encode algebraic properties of the tensors. Furthermore, these 
tensor networks can readily be defined on any lattice, allowing us to consider how 
the states represented by the tensor network vary with lattice geometry (in arbitrary 
spatial dimensionality) and the local on-site dimension. 

The states we consider display topological order. In Section 4 we show that they 
are the deconfincd phases of Abclian lattice gauge theories. The fractional quantum 
Hall effect [27] provided the first examples of topologically ordered phases. Topological 
phases of matter are also found in lattice gauge theories [28], spin liquids [29] and 
topological insulators [30] . Interest in these exotic phases of matter stems from their 
novelty of lying outside the conventional Landau symmetry-breaking paradigm. They 
are characterised by their long-range entanglement and topology-dependent ground 
state degeneracy. Their study is further motivated by the potential application of 
topological phases of lattice models to fault-tolerant quantum computers and quantum 
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memories. Kitaev's toric code [31] provided the first example of a topological quantum 
error-correcting code based on the Z 2 group. Our class of states is a subclass of those 
accessible within Levin and Wen's string- net formalism [32], and as such is already 
known to have tensor network representations which are fixed points of the TRG 
procedure [33, 34]. We provide a nascent and exemplary class of states, algebraically 
contractible via a considerably simpler route. 

2. Summary of results 

In the next section, we introduce tensor network states in a general setting, along 
with the special class of boolean tensors introduced in [35] and the pertinent rules 
regarding their manipulation. These are illustrated by the well known GHZ state. In 
Section 4 we briefly review models with topological order, then introduce Z 2 lattice 
gauge theory with a discussion of its physics as a case study. We refer to the "Z 2 
state" as the ground state of a Hamiltonian with Z 2 topological order. In Section 5 
we reproduce the tensor network state description of the Z 2 state from [33] (adapted 
for the hexagonal lattice) and present a new property of this TNS, namely that it is 
exactly contractible through the bialgcbra, Hopf and tensor fusion relations introduced 
in Section 3. We use our construction to calculate the two-point correlation functions 
of this state, and confirm that it has a zero correlation length. We also show a direct 
calculation of the topological entanglement entropy from the TNS. In Section 6 we 
discuss the excitations in Z 2 lattice gauge theory, and demonstrate how these can be 
included using "impurity" tensors in the network. The exact contractibility property 
enables the network to be contracted around the impurities provided their number is 
finite. We demonstrate this and consider string-like regions of the lattice containing a 
variable number of impurity tensors, calculating the reduced density matrix and von 
Neumann entropy for this region. As a second application, we calculate the topological 
entanglement entropy in the region around an impurity tensor. 

Section 7 takes the Z 2 TNS and generalises it to a class of tensor networks 
representing the deconfined phases of all finite Abelian lattice gauge theories. This 
is done in a way in which the tensors used to construct this class obey the same 
rules that enabled the exact contraction in the Z 2 case. We also show that these 
are, in a certain sense, the most efficient tensor network representations possible of 
these quantum states that can be expected. In Section 8 we translate these tensor 
networks onto square, kagome and triangular lattices and compare the properties of 
the corresponding quantum states. We find that on bipartite lattices the network will 
contract by precisely the same mechanism, with identical mathematical structures 
pointing to identical physics in every case. By contrast, on non-bipartite lattices the 
network remains contractible but with some important differences. In particular, these 
tensor networks are seen to represent non-trivial quantum states, but can only possess 
Z% topological order. Section 9 concludes the paper. 

3. Tensor network states and boolean algebra 

A tensor network provides an efficient means to describe a quantum many-body system 
of spins living on a lattice (or any system which may be mapped as such). Product 
states through to symmetry-breaking and topological phases are all accessible within 
this description. We use the physical lattice geometry to create a network in which 
tensors are associated to the physical spins. For each spin with coordination number 
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n, the corresponding tensor network description incorporates an order (n+1) tensor 
A™i at each vertex, with n internal indices a, /3, 7 . . . providing the network edges 
and one physical index r/tj representing the state of the spin in a physical basis. All 
internal indices are summed over in any network. Typically, internal indices are d- 
dimensional and the physical indices are D-dimensional, however in this paper we will 
always consider the specific case d = D. On a hexagonal lattice in which the spins are 
three-fold coordinated, a general tensor network state is expressed as 

l*>= E E A a^ A aL A Kr, A je\ • ■ • K m i»- W 

{ mi }=0 {a,/3,7...}=0 

The notation {m^}) indicates a quantum state of all spins on the lattice in which 
the i th physical spin takes the value m; in the chosen basis. Figure 1 gives a graphical 
representation of the quantum state in (1). In the figure, the square blocks represent 
tensors, open legs represent the physical indices, and lines connecting tensors denote 
the sums in (1) — contraction over pairs of indices. Periodic boundary conditions 
are typically employed on a sufficiently large network, so that the internal legs on the 
boundaries wrap around the lattice. In addition, requirements of translational and 
rotational invariance can often restrict the tensors to all be identical and rotationally 
symmetric, respectively, allowing for a highly compact state description. However, 
evaluating any probability amplitude or other property of such a state, such as reduced 
density matrices, involves contraction of a network. In general this is an unfeasible 
operation. To carry out the contraction naively, e.g. by contracting the tensors along 
lines in one direction first has exponentially large memory requirements. This is seen 
in extracting a single amplitude from a TNS on a square lattice: each tensor has d 4 D 
components and a line of n tensors contracted together forms a tensor with d 2n+2 D n 
components. This is prohibitive for all but the smallest of lattices. Numerical schemes 
to approximately perform this contraction exist (such as the MPS-MPO method [15] 
and TRG [16, 33]) but here our interest focusses on specific cases in which we can 
carry out the full contraction analytically for lattices of arbitrary size. 
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Figure 1. Graphical representation of a section of a general tensor network state 

with amplitudes f ({m;}) = Y,{aj,-y...}=0 A ™M A 2L A ™t n A ™L ■ ■ •< in which the 
physical degrees of freedom are arranged on the vertices of a hexagonal lattice. 
The square blocks represent tensors, the lines between them indicate summations 
over pairs of internal indices, and open (long-dashed) legs correspond to physical 
indices. 

We now briefly review the "constituent network components" introduced in [35] 
which will be the building blocks of our class of algebraically contractible states. We 
restrict the dimensionality of all tensor indices to d = 2 as relevant for spin-i systems, 
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and furthermore all tensor components to be or 1. This is the realm of the so-called 
"boolean tensors" [35, 36]. We introduce two tensors, COPY (•) and XOR (©), 

COPY XOR 



<5 'ijk 



1, i=j = k, 
0, otherwise, 



X 



ijk 



A 

1, i + j + k = (mod 2) , 
0, otherwise. 



(2) 



The legs of these graphical representations signify their indices. The COPY tensor 
defined above is a three-index generalisation of Kronccker's delta, and its extension 
to an arbitrary number of legs is straightforward. The second major component, the 
XOR tensor is related to the COPY tensor [37]. We can write it in terms of COPY 
through the contraction of each leg with a d = 2 Fourier matrix Hjk = (— l)- 7 ^, denoted 
by □, 




(3) 



We can likewise define n-legged XOR tensors in terms of n-legged COPY tensors. 
The XOR tensors have a number of useful diagrammatic properties. They are clearly 
fully symmetric under the exchange of legs. They have units, i.e. basis vectors which 
when contracted with any leg of the tensor give an identity matrix 




/A 



, (4) 



where |+) = (1 1) T and |0) = (1 0) T . These same vectors form copy-points with the 
two tensors 





(5) 



The COPY tensor obeys a so-called "tensor fusion rule"f allowing a network of COPY 
tensors to be amalgamated at will, e.g. 



(0) 





<J I 



f This is not related to the fusion rules in topological quantum field theory, which wc will refer to 
exclusively as "branching rules." 
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providing a purely diagrammatic interpretation for the equation Sijkl 
This is also true of the XOR tensor, J 





The COPY and XOR tensors together obey a bialgebra law [38, 39], 

m 



j 



V 



< 



JB- 



3^ 



(7) 



(8) 



Later on, we will consider a generalised tensor while maintaining the requirement 
that it obeys a bialgebra law with a higher dimensional COPY tensor. Expressing the 
bialgebra law (8) as a tensor equation, 

^ ^ Bjkn^lin ^ ^ 0~jpmO~krqBrnqiBp r l, (9) 
n mpqr 

we derive a condition that these generalised tensors Bijk must obey to satisfy a 
bialgebra law with COPY. All components Bjki = or 1, and further BjkiBjki = 
if i 7^ I. These place a heavy restriction on the number and location of the non-zero 
elements of B^k- Finally, the COPY and XOR tensors also obey a Hopf law [36, 39] 



-O- = H><F 



(10) 



We emphasise that these rewrite rules are defined up to an overall complex 
constant of proportionality, and that detailed proofs can be found in [35] . As a simple 
application of this graphical formalism, we can express the n-party GHZ state I^ghz) 
as 



I* 



GHZ 



)=2" 



(|00 . 

r 



.0) 



ii 



i)) 



(ii) 



It is straightforward to see that this represents I'I'ghz): contracting each physical leg 
with |0) or |1), we see that the coefficient is non-zero only when |00 . . . 0) and 1 1 1 ... 1 ) 
are inputted. We can calculate any two-site reduced density matrix by taking the 
complex conjugate of the state ((*ghz|) and contracting over all but two physical 
legs. In the chosen basis all amplitudes are real, and so the complex conjugate may 
be overlooked. Our convention is that upwards- and downwards-directed physical legs 
correspond to bras and kets respectively. The contraction of the state follows directly 
from the COPY tensor fusion rule. 



Pab = 



, — 1 — . 


» 1 


1 1 


' — T 


<, — 1 — . 


) 1 


1 ( 


. — 1 




J Later, we will consider a generalisation of this tensor to higher dimensions. Note that, unlike 
COPY, this generalised XOR only obeys a fusion rule if it is constructed from a direct sum of Z2 
groups. 
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-(|00)(00| + |11)(11|). (12) 

As noted above, the rewrites are only defined up to proportionality, so in (12) we 
have normalised at the final stage. Diagrammatically it is clear that p takes this form 
regardless of which sites we choose, showing that the correlation length for this state 
is infinite. 



4. Topological models 



Before introducing the tensor network state for the Z2 model and generalisations, we 
provide a brief introduction to the main aspects of Levin and Wen's class of string- 
nets [32], as our resulting quantum states will form a subclass of these models. A 
state in this class consists of a set of labels on the edges of a hexagonal lattice. The 
labels correspond to directed string types, and the Hilbert space is composed of all 
possible combinations of these labels on each edge. A particular string-net condensate 
is completely specified by a set of "data" [32] . 

(i) A list of directed string types, i = (0, . . . , N — 1). Each string type has a 
"conjugate" type — the same string with the reverse orientation, which is written 
as i* = j. Undirected strings have i* = i. 

(ii) A branching rule tensor = 1 if string types i, j and k may meet at a vertex, 
otherwise. We will use T to denote the branching rules tensor throughout, but 
in the Z2 case Ty* = Xijk- 

(iii) A quantum dimension di for each string type. The total quantum dimension T> 
is defined as V = (V.d?^ 1/a 



(iv) 



describing the recoupling relations between different 
string types. These are identities allowing for the amplitudes of locally different 
string-net configurations to be related to one another [32], viz. 



A six-index tensor Ff,j™ 



( 



V 




\ 



/ 



E 



V 




\ 



/ 



(13) 



In this notation, VP (Q) is the string- net amplitude for a configuration Q on the 
lattice. The piece displayed above is a local section of a string-net, the remainder 
is arbitrary and equal on both sides of the equation above. Fjfi™ satisfies the 
so-called "pentagon equation" to guarantee a single- valued wavefunction [40] . 



kp" n mns* Ikr* 



q* ki 



riq 
mis* 



(14) 



Many of the components of F^f™ are related by symmetry. The F-tensors contain the 



branching rules implicitly, for example 



T ijk - 



-F 



j*i*0> 



(15) 



Finding solutions to (14) is difficult in general, however it is the case that all groups 
and quantum groups generate a solution in which the string type index runs over the 
irreducible representations of the group and F^™ is its corresponding 6j-symbol [32] . 
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Figure 2. Graphical depiction of the tensor network for the Z2 state, from [33]. 
The network components are the COPY tensor S a p~f (•) and the XOR tensor 
X a pj (©)• The XOR tensors enforce the "closed-loop" constraint of the Z2 state. 

We now focus on Z2 lattice gauge theory as the simplest example of a topological 
quantum field theory. In terms of spin-i degrees of freedom living on the sites of a 
kagome lattice, its Hamiltonian is 

H = U T. (i-IrI-sE IK- j E^- ( 16 ) 

vertices v \ i£v / plaquettcs p i£p edges i 

The terms "vertices", "plaquettcs" and "edges" refer to those on the dual (hexagonal) 
lattice. We can also discuss this Hamiltonian using string language. The model is 
equivalent to one type of undirected string living on the edges of a hexagonal lattice. 
In terms of the eigenstates of er 2 , |0) and |1), we interpret the state |1) on an edge as 
having a string present, and |0) being empty. The strings correspond to the "electric 
flux" lines of the lattice gauge theory. Generally, we also assume g > 0, J > 0. 
This model is dual to an Ising model in (2+1) dimensions (in which the strings are 
domain boundaries), and contains a quantum phase transition [33]. When g = 0, 
the ground state is ferromagnetic with all sites spin-up. This is the confined phase 
of the lattice gauge theory, and in string language the ground state is absent of 
strings. As g/J increases past a critical value (g/J) c , the ground state becomes a 
weighted superposition of closed string configurations. In the limit J = 0, the above 
Hamiltonian (16) becomes exactly solvable, 

*=-£ ir-e ir ^ 

vertices v i£v plaquettcs p i£p 

Provided that the ratio g/U is finite, it does not affect the ground state, and so we 
set both U and g to 1 for simplicity. We focus on the ground state in this limit, the 
"Z2 state" which is an equal superposition of all configurations with closed strings, 

1*)= E 

S closed 

where <S* is a basis configuration of the lattice, and 'S closed' signifies the subset of such 
configurations in which there are only closed loops — no open or branching strings. 
This is the deconfincd phase of the lattice gauge theory. Despite its simplicity, this 
state is topologically ordered. As a final remark, the Z2 state is a stabiliser state, since 
it is an eigenstate of the set of stabilisers Yii^v °f with eigenvalue +1. 

5. Algebraic contraction of the Z2 state 

With the basic physics of Z2 lattice gauge theory now presented, we turn to a tensor 
network representation of an ideal Z2 state [33]. This tensor network is depicted in 
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Figure 2, and can now be seen to be composed of precisely the COPY and XOR 
tensors we have introduced in Section 3. Intuitively, the COPY tensors transmit the 
chosen basis states (|0), |1)) on the physical indices to the internal indices of the 
XOR tensors, which then enforce the branching rules (closed- loop constraint). Wc 
now demonstrate how the graphical rewrite rules introduced in Section 3 allow us to 
contract this state exactly. First consider an area in which all physical indices are 
contracted. The double-layered network structure results from contracting physical 
indices on two copies of the state (from (VP) and \^f)), 




Focussing on a small region of the lattice, the diagrammatic local rules previously 
introduced are put to use. We manipulate the network, using the fusion rule to 
rearrange the COPY tensors followed by the bialgebra law. 




= ^ (20) 

Consequently, the double layer network flattens to a single layer and can subsequently 
be contracted via the XOR fusion rule, 




We can use these rewrite rules to calculate quantities of interest, such as two- 
site reduced density matrices and the topological entanglement entropy. In a way 
analogous to that for the GHZ state, the method readily reveals the nature of the 
correlations in this state. Specifically the two-site reduced density matrix is given by 
a double-layer network similar to that above, but with the two sites of interest left 
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The local rewrites clearly expose the product state nature of the two-site reduced 
density operator. Consider the connected correlation functions (0(ri)0(r2)) c = 
(0(ri)0(r2)) — (0(ri))(0(r2)), in which r.i are positions on the lattice, and O is 
any local operator. Then (O(r)) = ±TrO and (0(n)0(r 2 )) = ±(TrO) 2 , hence 
(0(ri)0(r 2 )) c = 0. As such, the correlation length of this state is zero. This is 
also true of product states, yet this state is characterised by its pattern of long-range 
entanglement as revealed by the topological entanglement entropy [41]. We can employ 
a similar strategy to calculate this quantity directly. The entropy of entanglement S 
is defined as the von Neumann entropy of a reduced state pa of a bipartite system 
AB, 

S = -Tr (pa log pa), PA = Tr s (|*)(*|) . (23) 

Generally, the entropy of entanglement of a region in a topological phase scales with 
the region's boundary L as S = aL — 7 + . . ., with terms which vanish as L — > 00 
omitted [41, 42]. As the state has a correlation length of zero, we have precisely that 
S = aL—j. Here 7 is a subleading constant term which arises in the case of topological 
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phases (7 > 0) [41] . On a lattice, a region enclosing N sites has an ill-defined boundary 
L, so a scheme based on differences is employed to ensure the cancellation of the 
leading term. The Kitaev-Preskill topological entanglement entropy [41] is defined as 
Stop — $A + S B + S c - Sab - S B c - S C a + Sabc = -7, with A, B and C denoting 
mutually adjacent subrcgions of the lattice, such that the composite region A U B U C 
is simply connected as shown in figure 3. 



The zero correlation length of these states enables a simple calculation of the 
topological entanglement entropy, since the limit L — > 00 need not be taken. The 
leading term aL is cancelled in the construction, and all other potential L-dependent 
terms are zero. Consequently, we may proceed by considering just three adjacent sites. 
From the reduced density matrices above, Sa = log 2, Sab = 2 log 2. The three-site 
reduced density matrix pabc evaluates to 



This reduced density matrix does not factorise and is diagonal with 2 2 equal non- 
zero entries. It gives Sabc = 2 log 2 yielding the topological entanglement entropy 
Stop = — log 2 in complete agreement with [43] , but computed here by very different 
means. 

6. Excitations in Z2 lattice gauge theory 

The Z2 model has four kinds of elementary quasiparticle excitation operators: trivial, 
"magnetic flux", "electric charge" and "electric/magnetic bound states". The first three 
have bosonic statistics, while the last is fcrmionic. These operators act along strings 
on the lattice, creating pairs of quasiparticlcs, one at each end of the string. The 
excitations are gapped (with energy of order U, as in (16)). 

Within our formalism, an excitation may be included by inserting an "impurity 
tensor" ((§)) on a site which enforces an open string there. The rest of the lattice may 
be contracted around any impurities, leaving a network sufficiently small to contract 
directly. As an example, we calculate the density matrix of a string of n = 4 spins as 
shown in figure 4. The "open" tensors (O) are unspecified, and will be set shortly to 
be either standard Z2 or impurity tensors. 




Figure 3. Connectivity of regions as taken when calculating the topological 
entanglement entropy Stop = Sa + Sb + Sc - Sab - Sbc - Sca + Sabc- The 
regions and definition of Stop are chosen such that the term in S scaling with the 
boundary L cancels in Stop- 




(24) 
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Figure 4. Local perturbation to the ground state via the inclusion of impurity 
tensors. The hexagonal lattice is divided into a string-like region A, and the 
remainder region B. 



A contraction over the degrees of freedom in the bulk (region B) yields the density- 
matrix for region A: 



We can choose the impurity tensors to be intermediate between the extremes of open 
and closed loops, where 9 parametrises this set, 



Setting 9 = corresponds to the standard Z2 state, while 9 = tt/2 corresponds to 
enforcing open strings at the location of each impurity tensor. We can now set each of 
the (n + 1) tensors in (25) to be either the standard Z2 tensor or impurity X' tensor, 
and examine the variation of the von Neumann entropy with 9 and the number of 
impurities m. This is plotted in figure 5. 

Physically, if all of the vertices in the string are standard Z2 tensors, the reduced 
density matrix is proportional to the identity and so S takes its maximum value, 
S = n log 2. A second intuitive limit is to choose 9 = 7r/4 and set all (n + 1) tensors 
along the string to be the impurity X^ k . For this value of 9, the impurity tensor 
factorises, 



PA 




(25) 






(27) 



Substituting this into the density matrix (25) immediately yields 




* * 4 4 
^ $ ^ $ 

I + + + +)<+ + + + I 



(28) 
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SI log 2 



\ 
\ S 



m = 0, 1 



m = 2 



m = 3 



m = 4 



Figure 5. Von Neumann entropy S for the four spins in a string-like region of 
the lattice, in the configuration shown in figure 4. Each curve corresponds to 
varying the impurity tensor parameter 9 for a string with m of the five vertex 
tensors along the string set to the impurity Xi^, and the rest as the standard 
%2 Xijk. The impurity tensor X' is defined in (26) such that 9 = corresponds 
to the standard Z2 X-tensor, 9 = n/2 corresponds to enforcing open strings and 
9 = 7r/4 as the mid-point between these two extremes is an equal superposition of 
the two, and has a factorisation as in (27). Note, S depends only on the number 
of X ' tensors and the value of 0. It is independent of the relative positioning of 
these tensors. 



This four-spin subsystem is in a pure state — it is an island, unentangled with the 
remainder of the topological phase, a fact not easily revealed by other means. 

For general m, (n + 1 — m) log 2 < S < nlog2 for 1 < m < n + 1. For m = 
there are no impurity tensors and so no ^-dependence. The form of the curves can 
be understood as follows: if all tensors along the string are set to X' , 6 varies the 
subsystem between the standard Zi and a pure state. If instead n < m — 1 tensors 
are set to X' , provided that they are all adjacent it is possible to trace out the ends 
of the string leaving a shorter string entirely consisting of X' tensors. The degrees of 
freedom traced out were responsible for the entanglement and the residual non-zero 
S even at 6 — ir/4. If instead the X' tensors are not adjacent, this procedure cannot 
be carried out. Nevertheless (25) indicates that the moving the impurity tensor only 
rearranges the tensor components of p, and does so in such a way that S remains 
unchanged. 

As a second application, we can make use of the algebraic contraction property to 
study the topological entanglement entropy for the composite region ABC shown in 
figure 6 as a function of 9. This configuration consists of two impurity tensors placed 
arbitrarily on the lattice, and the regions A, B and C each containing one of the spins 
adjacent to the first impurity tensor. The overlap of this state with the "ideal" 
state is cos#. The tensor network for pabc contracts to 



PABC 




(29) 



= - diag (cos 6, sin" 9, cos 



, sin 
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The one- and two-site reduced density matrices arc proportional to the identity as 
before. Consequently, if we calculate the topological entanglement entropy as before, 
considering just these three sites, we find 

S top = - log 2 - cos 2 9 log (cos 2 6) - sin 2 9 log (sin 2 6) . (30) 

When we expand the region ABC to include more spins, we find (30) remains 
unchanged until the second excitation tensor is contained within the region ABC. 
At this point, we retrieve the result Stop = — log 2 as required for this phase. Thus, 
the excitation tensors introduce correlations between the local spins in the two regions, 
and this must be taken into account in the regions when calculating the topological 
entanglement entropy in this way. At both 9 = and 9 = ir/2 the three-site calculation 
gives the full value (—log 2) for the state; in the case of a 'pure' excitation, the no 
correlations are introduced because we have effectively just locally relabelled the spins. 
At the intermediate value 9 = 7r/4, the correlations are maximal. At this value of 9, 
X'ijk factorises to a product as in (27). 

If we place the second excitation tensor at the vertex immediately adjacent to 
the spin in region A (say), we have the situation considered in the first example and 
see that only Sa, Sab and Sac vary with 9. These do not cancel in Stop, leading to 
the 0-dependence in (30). If instead the second excitation tensor is not immediately 
adjacent to any of the spins in the region ABC, only Sabc varies, but precisely as 
required to yield the same result for Stop- 




Figure 6. Accessing excited states via the inclusion of impurity tensors. The two 
red starred vertices signify the impurity tensors, while the rest of the network is the 
standard Z2 state (figure 2). The partition of regions around one of the impurities 
is marked: the regions A, B and C contain one spin each, all immediately adjacent 
to an excitation tensor. 

Our formalism allows for straightforward computations concerning localised 
perturbations and excitations, and their effect on the von Neumann entropy of string- 
like regions as in figure 4 and the topological entropy as in figure 6. These are two 
examples of the physical quantities which are directly calculable by this method; the 
physics concerning such configurations of isolated excitations is not otherwise easily 
accessible. 

7. Generalisation to a class of contractible states 

We have shown that the Z2 state for a spin-i system has an efficient tensor 
network description which can be contracted efficiently by local rewrite rules to 
evaluate quantities of interest. We now generalise this network to provide an efficient 
description of all finite Abelian lattice gauge theories for spin-S 1 systems. We retain 
precisely the tensor network structure introduced in figure 2, but instead of using the 
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labels {0, 1} for the indices, we now replace them with the elements of the Z 2 group, 
{/, e} with e 2 = /. We can cquivalently define the COPY tensor as 

= { J ot^se!' («.A7)e{/,e} (31) 

In the case of the XOR tensor defined in (2), we make use of the group product. The 
branching rules T a ^ may be equivalently redefined as 

T ^-{ I otZrZL (°>A7)e{I,e} (32) 

We continue to represent T a ^ 7 graphically with ffi. The generalisation to all finite 
Abelian groups is straightforward, and proceeds by replacing Z 2 with a group G, for 
which (a, /3, 7) G G in the definitions above. We have restricted focus to Abelian 
groups, since in order to connect these to their string-net condensate descriptions, we 
will make use of the one-one correspondence between an Abelian group's elements and 
its irreducible representations [44]. In the case that G = Zd, Ta/37 is a very intuitive 
generalisation of XOR, namely a 'PLUS' tensor as introduced in [37]. A PLUS tensor 
with d-dimensional indices has coefficient 1 when its indices sum to (mod d) . In this 
case, it is connected to a d-dimcnsional COPY tensor through Fourier matrices as in 
(3), with components Hjk = e t2n: ' k ' d , where the integers j and k label the elements of 
the cyclic group. 




Diagrammatically, we represent H and its complex conjugate H* with Q and 
respectively. In the finite Abelian case, G can be decomposed as the direct sum 
G = Z ni ffi Z n2 © ... © Z„ k , in which n» are powers of primes. Equation (33) applies 
with H taking a tensor product structure H = H ni © H n2 © ... © H nk §. At this point, 
we have constructed a tensor network with an identical structure to before but with 
the tensors themselves generalised. When attached to a physical index, the COPY 
tensor can be thought of as converting the chosen basis states |z) to group elements 
on the internal indices of the T tensors. The COPY tensor sets the dimension of the 
physical degrees of freedom to be d = |G|, the order of the group, i.e. this is a spin 
S=±(\G\ - 1) state. 

We want these newly-defined tensors to be contractible by similar means to the 
Z 2 network. Our key requirements are (i) a connection of the T tensors to the COPY 
tensor through H, so that we can use the COPY tensor fusion rule to contract the 
bulk, and (ii) a bialgcbra law to flatten the lattice. The first of these has been 
demonstrated in (33), and the second holds as long as T a p 1 defined in (32) satisfy the 
condition (9) required for a tensor to obey the bialgcbra relation (8) with COPY. That 
this is true can be shown as follows. First, all elements T a p 1 = or 1. Second, the 
product of two elements T a j3^T a ps = 0, since T Q ^ 7 = 1 implies a(3j = I. The group 
structure guarantees the existence of unique inverses, so it follows that T , ct ^ 7 T Q ,^ < 5 = 0: 
if a/37 = It a fi& 7^ I unless 7 = 5. Then the condition is satisfied. 

Consequently, the contraction performed in equations (19-22) proceeds just as in 
the case d = 2, with only a slight alteration after the flattening step. In general, it is 



§ If instead we had based the tensor network around a non-Abelian group, this property would not 
hold and we expect the contraction of the network, if still possible, to be more involved. 
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no longer true that H is self-inverse. However, H in (33) can be freely replaced with 
H* and it is a property of Fourier matrices that HH* = I. Consequently, provided 
the lattice is bipartite, T a p 1 may be decomposed into COPY x H and COPY x H* 
on the two sublattices, respectively, 





(34) 



This enables the tensor fusion rule to contract the bulk state. We have therefore 
readily generalised the Zi state to a class of exactly contractible states based on finite 
Abelian groups. We can make some immediate observations about the properties of 
these quantum states. They are non-trivial as the T tensors exclude any configuration 
not satisfying the branching rules and further, they consist of an equal superposition of 
allowed configurations owing to the equality of all non-zero coefficients in the network. 

It just remains to show that these exactly contractible networks represent the 
deconfincd phases of Abelian lattice gauge theories. To do this we construct the string- 
net representation of these phases. The 6j-symbol ™ is defined only if the four triads 



of labels {i, j, to}, {k, I, to*}, {i, n, 1} and {j, fc, n*} are allowed by the branching rules. 

1 if 

otherwise. For 

completeness we note that these are the ground states of a very general class of exactly 
solvable Hamiltonians H [32], acting on d-dimensional spins on the edges of a lattice 



The 6j symbols of Abelian groups have all equal coefficients, therefore F^T 
the triplets above are allowed, and for convenience we set Fj^™ 



H 



E o.- E 



1 N 

a E B p- 



vertices v 



d 2 



(35) 



s=0 



plaqucttcs p 

Here Q v is a projector on a lattice site, projecting onto the subspace of allowed string 
configurations at that vertex, 

(36) 




and Bp acts on plaquettes p in a rather involved manner, 
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Note that Bp has a much simpler intuitive interpretation [32] , namely it acts on ^> to 
introduce a loop of type s around a plaquette p. This can then be incorporated into 
the lattice via the recoupling relations to give the form above. 

The branching rules can be extracted from the 6j-symbols via (15), and these are 
identical to the branching rules of our state because of the correspondence between the 
group's elements and its irreducible representations. The equality of the 6j-symbol 
components implies we can transform one string-net configuration into another with 
no alteration of the coefficient, thus the ground state consists of an equal superposition 
of all configurations consistent with the branching rules. This is precisely our tensor 
network state defined above. 

Finally for these states, we can calculate directly the topological entanglement 
entropy. The calculation proceeds as before, and again we find the three-site reduced 
density matrix pabc = J2 a p 7 T a /3j\a/3j)(al3j\, {a, /J, 7} = (l...d). From this and 
(32), the one- and two-site density matrices are again seen to be of product form. 
Combining these results, Stop = — logd for an Abelian lattice gauge theory with 
d = \G\. This is well-known: Abelian anyons have quantum dimension dk = 1, so for 
total quantum dimension T>, V 2 is just the number of superselection sectors = |G| 2 , 
e.g. 4 for Z 2 gauge theory. It is also known that the total quantum dimension relates 
to the topological entanglement entropy via Stop = — log 2? = — log \G\ [41]. However, 
it is demonstrated here using a very different method. 

We have taken the Zi tensor network and shown that it naturally generalises to a 
TNS class describing the deconfincd phases of Abelian lattice gauge theories. This TNS 
class also represents these states as efficiently as possible. A reduced density matrix 
pr> of a single spin having Hilbert space dimension D, has a von Neumann entropy 
which is bounded by < S < \ogD. Since an Abelian group of order \G\ generates 
a lattice gauge theory with a one-site reduced density matrix having S = log |G| , it 
is necessary that to reproduce the local physical properties of this state the tensors 
must have internal legs of dimension D = d > \G\. As our TNS class has internal 
indices with dimension d = \G\, it is therefore the most efficient representation of 
these quantum states in the sense that the internal bond dimension in the network is 
minimised. 



8. Dimensionality and coordination 

We have presented an efficient tensor network description of the deconfined phases of 
finite Abelian gauge theories on a hexagonal lattice. Through our construction, we can 
readily vary the lattice geometry on which the tensor network is based and determine 
the consequences for the underlying quantum state. No changes to the method are 
required on bipartite lattices, but for clarity we demonstrate contraction on the square 
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lattice. The overlap (^l^) exemplifies the contraction of the double- layer network, 




Focussing on a small region of the lattice, the diagrammatic rewrite rules are put to 
use again. For a tensor with components all in the set {0, 1} it is clear that (c.f. (20)) 




Formally, this follows from rearrangements through the bialgebra relation. The 
network flattens, and we can use the alternating decomposition into H and H* just 
as in (34) followed by the COPY fusion rule. 




That we can handle the tensor network here in exactly the same manner as on the 
hexagonal lattice demonstrates that this class of states possesses identical physics on 
these two lattices, and indeed on any other bipartite lattice such as the simple cubic 
lattice in three dimensions. In summary, this is well known, but demonstrated here in 
a particularly clean and simple way. 

We can also construct these tensor networks on non-bipartite lattices (e.g. kagome 
and triangular). We expect this to take us outside the string- net class, as the 
orientation convention for the branching rules [32] cannot be satisfied on all vertices 
of a non-bipartite lattice. The states remain exactly contractible, but demonstrating 
this requires more algebraic rewrites. On the kagome lattice, we first "flatten" the 
double-layer network just as before. Taking inspiration from an application of the 
TRG routine to treat the kagome lattice [17], we can rewrite the network as follows, 
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(41) 



Each PLUS tensor has been factored into a COPY tensor and Fourier matrices. We are 
free to insert two-legged COPY tensors, and combine these with the Fourier matrices to 
give two-legged PLUS tensors. We can then analyse the resulting triangular structures 
using group theoretic notions. Graphically we denote this structure as a white dot, 

\ a 



So 



lOt 



(42) 



To determine its properties, we input three group elements a,/3,7 € G into the right- 
hand side of (42). This tensor element S a/ 3 7 = 1 if a = , f3 = 7" 1 and 7 — a^ 1 
simultaneously. This implies that all indices must be equal, a = /3 = 7, and also 
that they must be self-inverse, a = or 1 . Otherwise S aj a 7 = 0. Hence this structure 
behaves as a COPY tensor but only for self-inverse group elements. We will refer to 
it as a "sub-COPY tensor." This tensor obeys a fusion rule both with itself and with 
the COPY tensor, 



(43) 



This allows for similar contractions to those previously considered. A particular 
difference arises in the topological entanglement entropy. For the kagome lattice, 
the four-site reduced density matrix is diagrammatically expressed as 






Pabcd — 




(44) 



This can be simplified further in special cases, in particular if G = (a direct sum 
of n Z2 groups) , then o = • and the T tensors can be combined via the fusion rule to 
give a four-site density matrix with the structure of (24). Otherwise, the effect of the 
sub-COPY tensor is that only self-inverse elements (equivalcntly, Z2 subgroups) make 
a contribution to Stop. By this, we mean that if G = Z2 © Z3 © Z4, we would calculate 
Stop = — log(2 x 2), as the Z2 group and Z2 subgroup in Z4 contribute, but the Z3 
group does not. This is entirely consistent with the expected complications arising 
from a simple translation of tensor network states for string-nets to non-bipartite 
lattices. The orientation convention for the branching rules cannot be satisfied, and 
this will have consequences for states containing directed strings. The exception to 
this is when G = Z^, in which case all strings are undirected and the full topological 
order is present. This does not preclude the possibility of constructing such theories 
on non-bipartite lattices, merely that the single-bond tensor network presented here is 
insufficient. Nevertheless, a state in our class on the kagome lattice constructed from 
the Z-i group (for definiteness) is in this case not topologically ordered (Stop = 0), yet 
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is an example of a non-trivial spin-1 quantum state with exact contraction properties 
as previously discussed. This can be compared with a state constructed from the Z 3 
group on a square lattice, which does yield a topologically ordered spin-1 quantum 
state. In both cases, each PLUS tensor enforces the same correlations in four nearby 
spins and it is then the network connectivity which sets whether or not this induces 
topological order in the corresponding quantum state. 

Turning finally to the triangular lattice, a similar result occurs. The network 
can be rearranged as follows, yielding a COPY and sub-COPY tensor network which 
contracts to the same result as on the kagome lattice. 




This is a fundamental difference between these tensor network states on bipartite 
and non-bipartite lattices. On bipartite lattices, we have efficient representations of 
finite Abelian lattice gauge theories. On non-bipartite lattices we still have highly non- 
trivial spin-S quantum states, but with reduced topological order. Another influence of 
the lattice geometry on the state, its coordination number, is revealed in the structure 
of reduced density matrices across different bipartite lattices. When we consider all 
the sites around a vertex, we see correlations in the reduced density matrix consistent 
with the branching rules. If we trace over just one degree of freedom, a product density 
matrix arises directly. 




(46) 



Consequently, on a bipartite lattice correlations can only be revealed in the reduced 
density matrix if all n spins around a vertex are considered. In the case of non-bipartite 
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lattices, tracing out one degree of freedom does not yield a a product density matrix, 
as the more complicated form of (44) suggests. 

9. Conclusion 

We have adapted well known laws from modern algebra and extended techniques to 
a known tensor network representation of the Z^ state, an archetypal topologically- 
ordered phase. In particular we have demonstrated the existence of simple graphical 
rules by which the network may be contracted, allowing for the straightforward 
computation of physical quantities of relevance. Via a natural generalisation in the 
on-site dimension of this network, we have found a class of tensor network states 
which display remarkably straightforward contraction properties as a consequence of 
the non-trivial algebraic properties of the tensor network components. This is to be 
contrasted against tensor network states in general, for which efficient contraction is at 
best approximate. Furthermore, we have shown that this class forms the most efficient 
TNS representation of finite Abelian lattice gauge theories. 

Our construction allows us to study the influence on the quantum state of varying 
the lattice geometry for a given tensor network. The identical mathematical structures 
on all bipartite lattices signify identical physics for each case. Furthermore, these 
tensor networks also yield non-trivial quantum states on non-bipartite lattices, but 
only the Zi subgroups of the base group G make a contribution to the topological 
entanglement entropy, in contrast to the bipartite case. The nature of the algebraic 
contraction is such that a finite number of impurity regions can be handled, as the rest 
of the network can be contracted around them, leaving a small network to consider 
directly. We have made use of this to calculate density matrices for string-like regions 
containing local perturbations, and shown how these can disentangle the region from 
the rest of the state. We have also used this to calculate how a perturbations induce 
correlations into the ideal Z2 state, which must be accounted for when calculating the 
topological entanglement entropy from the tensor network. 

There are numerous directions for expanding this work, with some of the more 
immediate generalisations including the study of states based on non- Abelian and 
continuous groups, and classical lattice models based on networks of this variety. More 
complex topological phases have tensor network representations, such as the double 
semion model and string- net states in general [33]. These tensor networks involve 
double- and triple-line structures, but may also have internal structure enabling a 
similar set of local algebraic rewrite relations. 
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